arXiv:1506.09046vl [math.AP] 30Jun2015 


The effect of a line with non-local diffusion on 
Fisher-KPP propagation 


Henri BERESTYCKI a , Anne-Charline C0UL0N b , 
Jean-Michel ROQUEJOFFRE^ , Luca ROSSI 0 

a Ecole des Hautes Etudes en Sciences Sociales 
CAMS, 190-198 avenue de France, F-75244 Paris cedex 13, France 
k Institut de Mathematiques de Toulouse, Universite Paul Sabatier 
118 route de Narbonne, F-31062 Toulouse Cedex 4, France 
c Dipartimento di Matematica, Universita degli Studi di Padova 
Via Trieste, 63 - 35121 Padova, Italy 


Abstract 

We propose here a new model of accelerating fronts, consisting of one equation with 
non-local diffusion on a line, coupled via the boundary condition with a reaction- 
diffusion equation in the upper half-plane. The underlying biological question is 
to understand how transportation networks may enhance biological invasions. We 
show that the line accelerates the propagation in the direction of the line and en¬ 
hances the overall propagation in the plane and that the propagation is directed by 
diffusion on the line, where it is exponentially fast in time. We also describe com¬ 
pletely the invasion in the upper half-plane. This work is a non-local version of the 
model introduced in [16], where the line had a strong but local diffusion described 
by the classical Laplace operator. 
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1 Introduction 

This work is a continuation of a program of investigation started in [16, 15, 17] aiming at 
understanding the influence of a line with its own diffusion embedded in a region where a 
classical reaction-diffusion process takes place. The underlying biological motivation is to 
understand how a line or a network of roads, lines and streams, may enhance biological 
invasions. Besides its obvious motivation from ecology for invading biological species, this 
type of questions, also arises in many other contexts, in cellular biology (for cell division 
mechanisms, see e.g. [30] for numerical methods , [5] for mathematical study by entropy 
methods), chemical engineering (crystal growth processes, see [ ]) or biophysics (see e.g. 

[ ] for a linear stability analysis) and medicine, sometimes in three dimensional regions 

bounded by a surface having certain properties different from the ones in the bulk. 
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There are many situations where networks of lines play a major role. In epidemiology, 
the spreading of certain diseases is known to be strongly dependent on the communication 
network. In the classical example of the “Black death” plague in the middle of the 14th 
century, the roads connecting trade fair cities allowed the epidemics to expand Northward 
at a fast pace. The account by [56] describes how it then further spread inwards from 
these roads to eventually cover whole territories. A recent invasive species in Europe, 
related to climate change, the Pine processionary moth has progressed at a faster rate 
than anticipated. An hypothesis in this context is the role played by roads in allowing 
jumps. See the update in the interdisciplinary volume [ ] on this issue and the related 

public health concerns. Another example of the effect of lines on propagation in open 
space comes from the influence of seismic lines on movements of wolves in the Western 
Canadian Forest. These are straight lines dug across territories by oil companies for 
the purpose of exploring and monitoring oil reservoirs. The study in [ ] reports the 

observation that populations of wolves tend to move and concentrate on seismic lines, 
allowing them to move along larger distances. 

In [16], [17], three of the authors introduced and studied a new model to describe 
invasions in the plane when a fast diffusion takes place on a straight line. In this model, 
the line {x = 0} in the plane M 2 - referred to as “the road” - carries a density u(x,t) 
of the population. The rest of the plane is called “the field”, and the density there is 
denoted by v(x,y,t). By symmetry, the problem may be restricted to the upper half¬ 
plane fi:=Kx(0, +oo), where the dynamics is assumed to be represented by a standard 
Fisher-KPP equation with diffusivity d. There is no reproduction on the road, where the 
diffusivity coefficient is another constant D. The road and the field exchange individuals. 
The flux condition results from the road yielding a proportion fi of u to the field, and a 
proportion u of v\ y _ Q jumping from the field on the road. The opposite of the flux for v 
appears as a source term in the equation for u. 

The system thus reads as follows: 

{ d t u — Dd xx u — uv\ =Q — /in, x e R, t > 0 
d t v - dAv = f(v), (x, y) G O, t > 0 (1.1) 

~dd y v\ y=0 = im- w\ y=Q , xel, t > 0, 

where d, D,n, u are positive constants and the function / is smooth and satisfies the 
Fisher-KPP condition: 

/(°) = /( 1 )=°> /> 0 in (0,1), / < 0 in (l,+oo), f{s) ^ /'(0)s for s > 0. 

The initial conditions are: 

?4 =0 = v o in u\ t=Q = u 0 in M. 

Let Ck denote the classical Fisher-KPP spreading velocity (or invasion speed) in the field: 

c K ■= 2\/ df'{{)). 

The fundamental paper [40] analyzed the solution to the problem 

Ut — d xx u = f(u), n(x,0) = H(x) (the Heaviside function). 


3 



Among other things, [40] shows that, modulo a shift in time, the solution converges to a 
travelling wave of speed Ck- This is also the asymptotic speed at which the population 
would spread in any direction in the open space - in the absence of the road - starting 
from a confined distribution, i.e. with compact support (see [2], [3]). 

The main result of [16] is the following 

Theorem 1.1 ([16]). There exists c*(D) ^ ck such that, if (v,u) is the solution of 
(1.1) emanating from an arbitrary nonnegative, compactly supported, initial condition 
(v 0 ,u 0 ) ^ (0,0), it holds that 

Vc > c*, lim sup \(v(x,y, t),u(x,t))\ — 0, 

t- S-+00 \ x \ >ct 

v>0 ( 1 . 2 ) 

Vc < c*, a > 0, lim sup \(v(x,y,t),u(x,t)) — {vfp, 1)| = 0. 

t^+OO \ x \ <c t 
0 ^y<a 

Moreover, c*(D) > Ck if and only if D > 2d, and 

C *{D) Cqq Vd as D —>■ +oo, with > 0. (1.3) 

In other words, the solution spreads at velocity c* in the direction of the road, and the 
propagation is strongly enhanced when D is large, even though there is no reproduction 
on the road. This theorem is completed in [17] by a precise study of the expansion in the 
held and the determination of the asymptotic speed of propagation in every direction. 

Theorem 1.2 ([ ]). There exists w* G C x ([—7 t/2, 7t/2]) such that 

Vy > u>*(V), lim v(x 0 + yisinV, y 0 + ytcosV,t) = 0, 

t— >+oc 

V0 ^ 7 < u>*(V), lim v(xo + ytsinV, y 0 + ytcosV, t) = 1, 

t —>+oo 

locally uniformly in (xo,2/o) € and uniformly on sets (y, V) G {M+ x [—tt/2, 7t/2], |y — 
tu*(i?)| > r]}, for any given rj > 0. 

Moreover, w * ^ ck and, if D > 2d, there is d 0 G (0,7t/2) such that w*(d) > ck if 
and only if |V| > Vo- 

This theorem provides the spreading velocity in every direction (sin V, cos d), and 
reveals a critical angle phenomenon: the road influences the propagation on the held not 
only in the horizontal direction, but rather up to an angle n/2 — d 0 from it. It is further 
shown in [17] that Vo —> 0 as D —» +oo. The theorem is illustrated by the numerical 
simulation of Figure 1, reported from the second author’s PhD thesis [25]. See also the 
review [ ] for a discussion of some of these aspects, as well as some open questions. 

Model (1.1) is an example of propagation guided by a set of lower dimension, a topic 
that has recently attracted much interest. See for instance [1] and [ ] for non-local 

models for front propagation guided by favored genetic traits, and [6] for the study of 
fronts guided by a line. In the case of interest here, the scenario depicted by Theorems 1.1- 
1.2 displays a propagation speed-up , or propagation enhancement: the spreading velocity 
can be much larger than the reference speed Ck, but remains asymptotically constant. 
One may wonder if different conditions may lead, not only to speed-up, but acceleration: 
in other words, the front velocity increases in time. In (1.1) we have only considered 
standard, local diffusion. So, it is natural to ask if non-local diffusion will accelerate the 
propagation, rather than only speed it up. 

The goal of this paper is to investigate this question. 
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Figure 1: Level sets of v 

2 Model with non-local dispersal on the road 

2.1 Description of the model and question 

To account for the possibility that individuals on the road may move much faster than in 
the held, we will consider a non-local diffusion on the road: this means that the underlying 
processes modelling the displacement of individuals are jump processes. When those are 
stable Levy processes, the corresponding diffusion operator is the fractional Laplacian: 
this is the choice we make here. Thus the system under study is 

{ d t v — Av = f(v), i6l, y > 0, f > 0 

d t u + (— d xx ) a u = —[M + vv — ku, i6l, y = 0, t> 0 (2.1) 

—d y v = yu — vv, i6l, y = 0, t > 0. 

We assume that the reaction term / is strictly concave and smooth, with /(0) = /(1) = 0. 

This assumption will always be understood in the following without further reference. 
Note that the second equation, representing the evolution of the density on the road, 
does not involve a reproduction term, but only possibly a mortality term — ku , with 
k ^ 0 constant. The exchange factor /i is a positive constant and we assume, without 
loss of generality, that v — 1. The operator (— d xx ) a , 0 < a < 1, is the fractional 
Laplacian 

f c\ \a / \ f U ( X ) ~ U (y), 

(~d xx ) “(a) : = c « [a ._ y |i +2 * d y- 
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The constant c a is adjusted so that {—d xx ) a is a pseudo-differential operator with symbol 
|£| 2a . Note that, when a < 1/2 this has to be taken in the sense of convergent integrals, 
while when a ^ 1/2 this definition should be understood in the principal values sense. 
We complete the system with bounded, continuous initial conditions n(-,-,0) = Uo and 
w(-, 0 ) = u 0 . 

The goal of the present paper is to understand how, and at what speed, the level sets 
of u and v will spread for large times. We will see that, in this system, propagation takes 
advantage of both the fast diffusion on the road and the mass creation term in the field. 


2.2 Main results and interpretation 

The first step is to identify a unique steady state of the system and to show that it is the 
global attractor of the evolution problem. 

Theorem 2.1. Problem (2.1) admits a unique positive, bounded stationary solution 
( V S ,U S ). Moreover, ( V S ,U S ) is x-independent, and solutions (v,u) to (2.1) starting from 
nonnegative, bounded initial data (uq,Mo) ^ (0,0) satisfy 


(v(x,y,t),u(x,t)) —> {V 8 {y),U a ), 

t— >+OG 


locally uniformly in ( x,y ) G M x [0, +oo). 

This Liouville-type result is a consequence of [15] (see also [ ]). Indeed, it is shown 

in Lemma 2.3 there that the solution (v,u) to (1.1) lies asymptotically between two 
positive x-independent stationary solutions (Vi ,Uf) and (V 2 , ) - One can check that 
this property is proved without exploiting the equation on the road, and hence it holds 
true for (2.1). Notice that (V\,U\) and (V 2 , £^ 2 ) are solutions to (2.1) as well, because 
they are x-independent. Then, Proposition 3.1 of [ | implies that (2.1) admits a unique 

positive x-independent stationary solution (V s , U s ), concluding the proof of Theorem 2.1. 

In the case without mortality ( k = 0) we have the trivial solution {V s , U s ) = (1, 1/p) 
(recall that v = 1). In general, V s = V s (y) and we know that V s (+ 00 ) = 1 . 

The issue is now to track the invasion front, and this is done in the next two theorems. 


Theorem 2.2. (Propagation on the road). Let (v,u) be the solution of (2.1) starting 
from a nonnegative, compactly supported initial condition (uq,Mo) ^ (0,0). Then, setting 


7 * : = 


m 

l + 2a’ 


1. V 7 > 7 *, lim sup u(x,t ) = 0, 

t- 5- + 0 O 13 .| ^ e7 t 

2. V 7 < 7 *, lim sup | u(x,t) — U s \ = 0. 

t- s>+oo 

Thus, the process at work on the road is the same as that of an effective reaction- 
diffusion equation of the form: 


d t u + ( -d xx ) a u = f e ff(u), ( 2 . 2 ) 

where /e//(0) = f'( 0). In particular, spreading will occur at the same rate as was com¬ 
puted in [23 1 for this equation (see also [ 22 ]), despite the fact that there is no reproduction 
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on the road. Thus, the road has once again a dramatic effect, and we have therefore iden¬ 
tified a new mechanism for front acceleration. 

Turn to the propagation in the field. Recall that Ck = 2a//'( 0) is the KPP velocity. 

Theorem 2.3. (Propagation in the field). Under the assumptions of Theorem 2.2, for 
all 6 e (0,7r), we have: 

1. \/c> Ck/ sin0, lim supu(rcos^-rsin^t) = 0, 

2. Vc < cx/sind, lim sup \v(r cos 9,r sin 9, t) — V s (r sin 6)\ = 0. 

t-f+oc o ^.r^ct 

The speed of propagation in the direction (cos(0), sin(0)) is thus asymptotically equal 
to Ck/ sin(0). When 6 is close to 0, this speed tends to infinity, which is consistent with 
Theorem 2.2. And so, the front is, asymptotically, close to a straight line parallel to the 
x-axis moving vertically at velocity Ck- One may interpret it as follows: the invasion on 
the road is so fast that the whole system behaves just as if the density in the field only 
saw the condition v = 1 at the boundary, as in the effective equation Vt~ Av = f(v), with 
v(x, 0,t) = 1. It is easy to see that this gives the correct behavior, by trapping v(x,y,t ) 
between two suitable translates of the solutions of 

d t v - d yy v = f(v) for y e K, v(y, 0) = (1 ± e)H(y), 

for every small arbitrary e. The whole scenario is summarized by the numerical simulation 
of Figure 2, once again reported from [25]: 

2.3 Organization of the paper 

First, we give a brief overview of the context of the problem we study here in Section 3. 
There we present a review of the existing literature related to our model: propagation 
enhancement, propagation acceleration, as well as further results obtained on the model 

(1.1) . The proof of the exponential in time asymptotic spreading in our model starts 
in Section 4, where we outline the strategy of the proof of Theorem 2.2. The idea is 
to trap the solution between a supersolution and a subsolution whose level sets move 
asymptotically at the same speed. This turns out to be a rather delicate task, the 
construction of the supersolution being quite different from that of the sub-solution. The 
former is achieved in Section 5, together with the main computations. Section 6 is devoted 
to the construction of an auxiliary subsolution for a ID transport equation. This is the 
building block used in the rather computational Section 7 to obtain the subsolution to 
the full system. Some weaker versions of the second statement of Theorems 2.2 and 2.3 
- with the convergence to U s and V s replaced by a positive lower bound - are derived 
in Sections 7.4 and 8 respectively. These bounds are used in Section 9 to complete the 
proofs of Theorems 2.2 and 2.3. Finally, an appendix studies the Cauchy problem for 

(2.1) and provides regularity and a comparison principle, similar to that of [16]. 

3 A review of front speed-up and acceleration 

In this section, we present a general overview on front propagation enhancement, which 
has seen an important development in the last 15 years. We start with the case of a single 
local equation. In the second subsection, we present further results obtained on (1.1). In 
the last one, we give an overview of the mathematical literature on accelerating fronts. 
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Figure 2: Level sets of v 

3.1 Front propagation and speed-up: an overview 

It is well known that diffusion, when coupled with reaction, gives rise to propagating 
fronts. The most common situation is the development of fronts travelling at constant 
speed. The basic result concerns front spreading in a homogeneous medium. Let us recall 
its main features. The equation reads 

u t — A u = f(u), u | ^ 0, ^ 0, compactly supported. (3.1) 

We are looking for a function R : M + — > M + such that 

Ve > 0, lim inf u(x.t) > 0, and lim sup u(x,t) < 1. (3.2) 

t-H-oo |*[<R((l- e )t) *-»•+«> |*|>«((l+€)t) 

A fundamental result in this direction is that of Aronson-Weinberger [3], which accounts 
for how a front develops from a compactly supported datum. 

Theorem 3.1 ([ ]). Letu(x,t) be the solution of (3.1). Then, still setting ck '■= 2a//'(0), 

1. Vc > ex, lim sup u(x,t ) = 0, 

t_s> +°° \x\^ct 

2. Vc < ck, lim inf u(x,t ) = 1. 

t—l-oo |a:|<ct 
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Figure 3: Transition between 0 and 1 

In other words, (3.2) holds with R(t) = cr-£. An important literature has developed 
from there; to discuss it is out of the scope of this paper. Let us just mention the 
forthcoming book [8], which provides a complete overview of the question. 

In the presence of heterogeneities, quantifying propagation can be quite difficult. 
The pioneering work in this held goes back to the probabilistic arguments of Freidlin 
and Gartner [ |. They studied KPP-type propagation in a periodic environment and 

showed that the speed of propagation is no longer isotropic: propagation in any direction 
is influenced by all the other directions in the environment. They provide an explicit 
formula for the computation of the propagation speed. Many works have followed since 
then, we mention for instance the references [10] for a definition and estimates of the 
spreading speed in periodic environment, [13] for the particular case of space dimension 1, 
and [9] for a general definition in arbitrary heterogeneous media, as well as new estimates 
of the heat kernel. 

Reaction-diffusion equations in heterogeneous media since then is an active held and 
the question how the environment may enhance propagation has received much attention. 
The hrst paper in the domain is that of Audoly, Berestycki and Pomeau [ ], which studies 
models of the form 

d t v + AV ■ ( vq(x )) = Av + f(v), t e M, (x, y) G M x R N (3.3) 

where A > 1 is large and q is an imposed divergence free how held. They propose, by a 
formal analysis, various asymptotics for the propagation velocity, in particular if q(x) is a 
shear how q(x) = (a(y), 0) or a periodic cellular how held. This has triggered a series of 
mathematical works; let us hrst mention a general deep estimate by Constantin, Kiselev, 
Oberman and Ryzhik [24] of the bulk burning rate of (3.3), defined as: 

V(t)= I f(v(t, x))dx. 

Jr n 

This quantity turns out to be a relevant one, especially in problems where the how is 
time-dependent or strongly heterogeneous - thus precluding the reduction of the problem 
to travelling or pulsating waves. The study of speed-up properties of propagation by 
an advecting velocity held is continued, in a mathematically rigorous fashion, in [39]. 
This work derives upper and lower bounds for (3.3) in terms of A: V{t) ^ CA for shear 
hows, and V(t) ^ CA l A for cellular hows. When the how is a steady shear how, a 
recent paper of Hamel and Zlatos [36] makes the Kiselev-Ryzhik upper bound sharp: 

under a Hormander type condition on a, there exists 7*(a, /) ^ / a(y)dy such that 

Jt n ~ 1 
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the velocity c*(Aa,f) of travelling fronts of (3.3) satisfies lim 

A-»+oo 


c(Aa, f) 
A 


= 7 * 0 ,/)- 


Moreover the quantity 7* is the unique admissible velocity for the following degenerate 
system with unknowns 7 and U: 


A y U + (7 - a(y))d x U + f(U) =0(lx T w_1 ) 


lim U(x,y) = 0, 

x— H-oo 

lim U(x,y) = 1. 

X—> — OO 


(3.4) 


A similar situation will be discussed in the next section. 

For reaction terms / which are of the Fisher-KPP type, there is a relation between 
the principal eigenvalue and the amplitude of the velocity field, which has motivated the 
paper by Berestycki, Hamel and Nadirashvili [9] on principal eigenvalue problems of the 
form 

—A (p + Aq(x).V(j) = y(A)(j>, 

y(A) being the eigenvalue that is sought. The link between the existence or nonexistence 
of a first integral for the flow, and the size of fi(A), is elucidated. 


3.2 Speed-up by the line of fast diffusion: recent results 

The qualitative properties of model (1.1) have been studied further, with the discovery 
of new effects. In what follows, we list some natural questions and the answers that have 
been given. 


3.2.a Is the bound c*(D) ~ CooVD robust? 

Whether this asymptotic persists if the nonlinear term / is changed to a different type of 
source term is not obvious. Indeed, a priori trivial question; this could indeed be thought 
of as a special property of models with Fisher-KPP type nonlinearities. Moreover, the 
asymptotic for the velocity is obtained through algebraic computations. So, it is quite 
natural to ask whether this result persists with more general nonlinearities, where explicit 
computations are not longer possible. And, indeed, the property seems to be general, 
as shown by Dietrich [28], [27]. Consider the situation where fl is the strip M x (0,L), 
with Neumann boundary condition at y = L. The function / vanishes on an interval 
[0, 9\, and is positive on (0,1). Then (1.1) admits a unique travelling wave velocity c*(D). 
Moreover, the velocity still grows like y/D: we have c*(D) ~ where > 0 is the 

unique c such that the following problem admits solution: 


cd x 4) - d xx <i> — if) • fief), 
cdjp - ddyy'i/j = /(VO, 

< —ddyip — luj) — ip, 

dyi> = 0 , 

Jtf,<f>)(-°°,y) = (0,0), (V,0)(+°o,y) 


x E R, y = 0 

(x,y) e ft 

iGK, y — 0 (3.5) 

iGK, y = L 
(l,l//x), ye[0,L\. 


This is reminiscent of the Hamel-Zlatos situation, although the speed-up mechanism is 
quite different. 
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3.2.b How is the spreading velocity modified if additional effects are in¬ 
cluded? 

Another natural question is what happens if nonuniform transport effects are included. 
In the following model, we assume that on the road there is a constant transport q as 
well as a constant mortality rate k. The conditions in the field remain unchanged. The 
equations for u and v thus read: 


{ d t u — Dd xx u + qd x u — uv — yu — ku x G M, y = 0, t > 0 
d t v — dAv = f(y) (x, y) G O, t > 0 (3.6) 

— ddyV = jin — vv x e M, y — 0, t > 0. 


In [15], it is proved that Problem (3.6) admits asymptotic speeds of spreading cf in the 

D k q 


directions ±e\. Moreover, if — ^ 2 + 

d /'(0) 


T 


then cf = Ck, else cf > Ck- 


For definiteness, let us 


v'Wm’ 

Let us give a brief biological interpretation of this result, 
consider propagation to the right (that is, in the direction ei). 

First we see some expected effects: mortality on the road makes speed-up more dif¬ 
ficult by raising the threshold for D/d past which the effect of the road is felt. On the 
contrary, a transport q > 0 on the road facilitates speed-up (to the right) by lowering the 
threshold. When q < 0, the threshold is raised by the same factor. 

A less expected consequence is the following. In the absence of mortality on the road, 
the threshold condition reads 


D 

~d 


> 2 - 


q 

vW 



i.e., 


q> c K 



In other words, a transport q larger than ck speeds up propagation, no matter what the 
diffusivity ratio is. Biological situations are, for instance, the spreading of a parasite by 
a river, as reported for instance in [38]. For values of q less than c^(l — D/2d) - and, in 
particular, for large values of — q - spreading towards right occurs at the KPP velocity: in 
biological terms, propagation upstream against a river flow remains unaffected, whereas 
downstream propagation, in the direction of the flow, can be strongly enhanced. This 
yields a remarkable assymetry which, as a matter of fact, is also felt in the field: the 
asymptotic shape of the front may deviate significantly from that computed in the absence 
of flow field (c.f. [17]). 


3.2.c Nonlocal exchanges between the field and the road 

Up to now, we considered exchanges between the road and the field taking place in 
the infinitesimal vicinity of the road. It is therefore a natural question to examine the 
influence of the range of exchanges between the road and the field, in other words what 
happens when the exchanges are nonlocal. This was recently investigated by Pauthier, 
[51], [49] and [50]. The model under study is 


d t u - Dd xx u = / u(y)v(x, y, t)dy - u 

Jr 

d t v - dAv = f(v ) - v(y)v + y(y)u, 


11 


y(y)dy, i6l, t > 0 

(x, y) G M 2 , t > 0 


(3.7) 







The functions /i and v are smooth, L 1 functions: individuals can jump away from the 
road, but only very few of them can go very far from it. Notice that the initial model 
(1.1) is retrieved from (3.7), at least formally, by setting 


H(y) \= n5 y= o, u(y) = v5 y=0 


where Ji and V are two positive given constants. The question is how the spreading 
velocity is modified by the introduction of this non locality. Pauthier shows that the 
thresholds are quite stable, but also discovers surprising effects. 

1. Propagation enhancement. The threshold D = 2d is still there [ ]: if D ^ 2d, we 

have c*(D) = Ck; if D > 2d, then c*(D) > Ck■ Moreover, there is Coo > 0 such 
that c*(D) ~ Coqa/ZI as D —» +oo. The dynamics when // and v are close to Dirac 
masses is also quite stable [ ]: if (// £ ,z / e ) e>0 converge to (p5 y=0 ,v5y =0 ), then (i) 

cl(D) —* c*(D) as e —* 0 , and (ii) the limits t —* +oo and e —>• 0 commute. 

2. Variation of the spreading velocity with the range of // and v. How does c*(D) then 
vary when /i and v vary, their masses being respectively kept equal to Ji and Vl (i) 

1 y y 

a new threshold ([ ]). Take (fJ>(y), v{y)) = — (/^o(—), /^o(— )), the functions /i 0 and 

R R R 

uq having masses Ji and V. Then, if c*(D, R ) is the spreading speed, it holds that 



lim cJD,R)>ck 

R^+oo 


(ii) Moreover (see [51]), contrary to what intuition suggests, the spreading velocity 
is not always maximized - under the constraint that the masses of /i and v are kept 
equal to the constants ji and V - when exchanges are localized on the road. 

3.2.d Other effects 

The speed-up effects observed in the original model (1.1) are displayed in other various 
situations. Consider, for instance, a strip bounded by the road on one side and with 
Dirichlet boundary conditions on the other. In this case, Tellini [58] proves the existence 
of an asymptotic speed of propagation which is greater than that of the case without 
road and studies its behavior in the limits D —> 0 and D —y +oo. When the width of 
the strip goes to infinity, the asymptotic speed of propagation approaches the one of the 
half-plane model ( 1 . 1 ). 

Equations (1.1) have also an interest in higher dimensions, where they arise as models 
in medicine. One motivation would be to model the diffusion of a drug within a body 
through the blood network. The analysis of an TV-dimensional model is achieved by Rossi, 
Tellini and Valdinoci [55]. The authors consider a circular cylinder with fast diffusion at 
the boundary, which reduces to a strip between two parallel roads in the bidimensional 
case. The picture obtained in [ ] is similar to the ones described before: enhancement 

of the spreading velocity occurs if and only if the ratio between the diffusivities on the 
boundary and inside the cylinder is above a certain threshold. The authors investigate 
the dependence of the spreading velocity with respect to the radius R of the cylinder, 
discovering that it is monotone increasing if the ratio between the diffusivities is below 
2 , whereas, if the ratio is larger than 2 , the dependence is no longer monotone and there 
exists a critical radius R = Rm maximizing the velocity. 
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Let us end this review with the recent preprint [34] by Giletti, Monsaingeon and Zhou, 
which extends Theorem 1.1 to the situation where /j and v are replaced by 1-periodic 
functions. 


3.3 Front acceleration: a review 

As opposed to the previous situation, non-local diffusion may cause acceleration of fronts. 
This phenomenon has also long been identified. In the context of ecology, Kot, Lewis and 
Van den Driessche [42] study, both numerically and heuristically, discrete time models of 
the form 

/ +oo 

k(x - y)f{N t {y))dy = (k * N t )(x). (3.8) 

-OO 

When the decay of the convolution kernel k is slow enough, the authors observe accel¬ 
erating profiles rather than travelling waves. Similar properties have been noticed, still 
from the numerical point of view or in the formal style, for models with continuous time, 
e.g. reaction-diffusion equations of the form: 

u t + Lu — f(u), t > 0, x £ R n (3.9) 


where / is of the Fisher-KPP type, and L a non-local diffusive operator. Typical examples 
are Lu = (—A)", or Lu = k * u — u (notice that, with this last kernel, (3.8) is the exact 
analogue of (3.9)). See [46] for a rather complete review. The basic heuristic argument 
for acceleration is the following: since / is concave, a good approximation of the dynamics 
of (3.9) is given by that of the linearized equation; this entails studying the level sets of 

v(x,t) = e^'^ t e~ tL u Q (x). 


In the case L = (—A) Q , and for a compactly supported datum u 0 (x), we have e tL ~ 
-— , yielding that the level sets of v spread like e^'^ >t ^ N+2a ^ . More generally, if e~ tL 

CC ^aol 

decays spatially slower than any exponential, this yields accelerating level sets. 

Mathematically rigorous proofs of acceleration, and precise identifications of the mech¬ 
anisms responsible for acceleration, are more recent. The first paper in this direction is 
that of Cabre and the third author [ ] for L = (—A)“, which proves that the level sets 

of (3.9) spread asymptotically like e-f , ( 0 ) i /( jV + 2o A 


Theorem 3.2 ([ ]). Let u(x,t) be the solution of (3.9) with L = (—A)“ ; starting from 
a compactly supported initial datum uq ^ 0, ^ 0. Then we have: 


1. Vc > 


m 

N + 2a 


lim sup u(x, t ) = 0, 

t-H-oo | x \^ e ct 


2. Vc < 


m 

N + 2a 


lim inf u(x,t ) 

£—>■+00 |o;|<e c£ 


1 . 


Here, (3.2) holds with R(t) = fL'^ t ^ N+2a ' > . See [ ] when e~ tL is a general Feller 

semigroup. A first question of interest is what happens as a —> 1, or how to reconcile 
Theorems 3.1 and 3.2 in the limit a —> 1. The second and third author address this 
question in [26]: for a close to 1, propagation at velocity 2+ fiff) occurs for a time of 
the order |Ln(l — a)|; from that time on, Theorem 3.2 applies. 
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When L is of the convolution type, the relevant result is that of Gamier [32], who 
proves: (i) super-linear spreading for (3.9) as soon as the convolution kernel k decays 
more slowly than any exponential, (ii) exponential spreading when k decays algebraically 
at infinity. The precise exponents are not, however, given. Accelerating fronts can be 
observed in (3.9) even when L = —A: Hamel and Roques prove, in [ ], that it is enough 

to replace the compactly supported initial datum with a slowly decaying one. This paper 
is also the first to identify, in an explicit way, that the correct dynamics of the level sets 
is given by that of the level sets of the ODE 

ii = f{u), u(0,x) = u 0 (x). 

Whether or not Theorem 3.2 is sharp is a natural question. One may indeed wonder 
whether the exponentials should be corrected by sub-exponential factors. Cabre, and the 
second and third authors prove in [21] that the exponentials are indeed sharp, in other 
words that any level set is trapped in an annulus whose inner and outer radii are constant 
multiples of (A (o )t/(N+2a)_ p or that, they devise a new methodology which extends to the 
treatment of the models of the form 

u t + (—A ) a u = fi{x)u — u 2 , (3.10) 

with p > 0 and 1-periodic. Surprisingly enough, the invasion property (3.2) can still 
be described by the single function R{t) = e A ° t /( N + 2a ) j where (—Ao) is the first periodic 
eigenvalue of (—A)" — p(x). The method consists in two steps: (i) one shows that w(-, 1) 
decays at the same rate as the fractional heat kernel, (ii) one constructs a pair of sub 
and supersolutions that have exactly the right growth for large times and the right decay 
for large x. This proved to be a more precise approach than all the previous studies, 
which mainly relied on the analysis of the linear equation. This mechanism, which is 
quite different from what happens in the case a = 1, was later on described by Meleard 
and Mirrahimi [48], with a different viewpoint. Their work is in the spirit of the Evans- 
Souganidis approach for front propagation [31]; to take into account the fact that the 
propagation is exponential, the authors modify the classical scaling ( x,t ) i-> (t/e,x/e) 
into (x,t) (t/e, l^l 1 / 6 ); they apply the Hopf-Cole transformation to the new equation 

and derive a propagation law for the level sets. 

The analysis of [21] can be pushed further, to prove that in fact a strong symmetriza- 
tion phenomenon is at work, see [ |. The result is the following: when uq is compactly 

supported, then the level sets of the solution u(x,t ) of (3.9), with L = (—A)", are 
asymptotically trapped in annuli of the form 

{g[u 0 ]e /,(0)t/(JV+2a) (l - Ce~ st ) < |x| ^ q[u 0 ]e f '^ N+2a \ 1 + Ce~ 5t )}. 

Here, q[u 0 ] > 0, C > 0 and 5 G (0, f(Q)/(N + 2a)) are constants depending on u 0 . 

Let us end this review by mentioning a different mechanism of acceleration in kinetic 
equations. Here, an unbounded variable is responsible for acceleration of the overall 
propagation. A first model is motivated by the mathematical description of the invasion 
of cane toads in Australia, ft has the form 

d t n — a(9)A x n — DA g n = n( 1 — p n ), t > 0, x G M. N , 9 G 0 (3.11) 

where n{t, x, 9) is the density of individuals, and 9 a genetic trait. The quantity p n (x, t ) is 
the integral of n over 0. The coefficient a(9) >0 may be unbounded, as well as the state 
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space 0. This influences the dynamics of (3.11): when 0 is unbounded and a(9) = 6 , 
the note [19] gives a formal proof that the level sets of n develop like f 3//2 . This, by the 
way, can also be seen by computing the fundamental solution of the of (3.11), linearized 
at n = 0, and is related to the previous heuristics concerning the fractional Laplacian. A 
work under completion by Berestycki, Mouliot and Raoul [18], gives a rigorous proof of 
the computations of [19]. 

A related model is the BGK-like equation 

d t g + v ■ V x g = (M(v)p g - g) + p g (M(v) -g), t> 0,xeR,veV (3.12) 

analyzed by Bouin, Calvez and Nadin in [20]. The underlying biological situation is that 
of a colony of bacteria. The unknown g(t, x,v ) is the density of individuals, and v a 
velocity parameter. The set V may be unbounded, and this, as before, influences the 
dynamics of (3.12). The function M(v) is a reference velocity distribution. The quantity 
p g (x,t ) is the integral of g over V. The Fisher-KPP equation (3.1) arises as a limiting 
case of (3.12) under the scaling (t,x,v) i—> (t/e 2 , x/e, £ 2 v), as e —* 0. When M(v) is a 
Gaussian distribution, a level set of a solution g(t,x,v ) starting from an initial datum 
of the form g(0,x,v) = M(v)ti x<XL y is trapped in an interval of the form(cit 3 / 2 , C 2 ^ 2 ) 
where q are universal constants. This result is obtained by the construction of a pair of 
sub/super solutions of a new type. The authors conjecture that acceleration will always 
occur when V is unbounded. 


4 Strategy of the proof of Theorem 2.2 and comments 


4.1 The main lines of the proof 

A first idea would be to try to adapt a new, and quite flexible, argument devised by the 
second and third authors of the present paper, together with X. Cabre, in [21]. For the 
equation 

u t + (—A ) a u = f(u), t > 0, x G 
0 < a < 1 , with u( 0 ,.) compactly supported, they prove the 


Theorem 4.1 ([21]). We have, for a universal C > 0: 

C~ x , , C 


1 + e 


—Kt \x\ N+2a 


^ u(x, t) ^ 


1 + e 


-ret \x\ N+2a 


(4.1) 


This is done by introduction of the invariant coordinates £ = xe xt , and a pair of 
sub/super solutions is sought for in those coordinates, where u solves 


d t u - A£ ■ V 5 u + e~ aXt {-A) Q u - f[u) = 0 


(4.2) 


The construction of sub/super solutions for (4.2) is connected to the existence of suitably 
decaying solutions of 

on the whole line, which just amounts to finding entire solutions of if) — /('if). Try¬ 
ing to extend this idea here amounts to rescaling the x variable, defining the functions 
v(£,y,t) := v(e' yt f, y, t) and «(£, t) := w(e 7t £, t), with the idea that 7 = /'( 0 )/(l + 2a). If 
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we - formally - neglect the diffusive terms e 2 lt v^ and e 2 ait (—d^) a u, that should go to 
0 as t tends to +oo, we end up with the following transport system 

f d t v-jid^v-dyyv = f(v), £ G R,y > 0,t > 0, 

< d t u —'yffd^u = —/. m + v — ku , £ G M, y — 0, t > 0, (4.3) 

[ — d y v = im — v, (6l,|/ = 0,t>0. 

The idea is to look for stationary solutions to that system, and try to deform them. 
However, we are not able to carry out that program, and there is a deep reason for that. 
The subsolution will be constructed in a different way than the supersolution, which will 
result in a loss of precision in estimating the propagation speed on the road. 

1. The upper bound. We use the classical remark that f(v) ^ /'( 0)u to bound the 
solution of ( 1 . 1 ) by that of the linearized system at v = 0 , i.e. the solution (v,u) of 


{ d t v — Av 
d t u + ( -d xx ) a u 
-d y v 


f(0)v, i6l,|/>0,i>0 

—fm + v — ku, i6R,|/ = 0,i>0 
/in — v, x G 1, y = 0, t > 0, 


What will be of interest to us will be the behavior f u on the road, the rest of the 
solution being handled with standard arguments of parabolic equations. The main 
results that we will prove is the existence of a constant c a > 0 such that 


u(x,t ) 


C(y& 




(k + / , (0)) 3 |x| 1+2 "f2 


as 


|a:|, t —y +oo. 


(4.4) 


This will give, for all 7 > 7 * = 


m 

1 + 2 a 


lim u(x,t) = 0 uniformly in |a;| ^ e 7i . 

t —^-)-oo 


In fact, we have, for t large enough 


/'( 0 ) 


{x e M j u(x,t ) = A} C < x e M | |x| ^ C\t 2 (1 + 2“) e(i+ 2o! ) 


2. Lower bound. We apply the methodology introduced in [21], but we adapt it in 
an important fashion: since it seems difficult to construct a stationary subsolution 
to the rescaled transport problem, we work in a strip of width L instead of the half 
plane and let L go to infinity. An explicit subsolution is constructed under the form 


v(x,y,t) 



if 0 < y < L(1 — |) 


u(x,t) = Ch4>{xBe 7 *) , 


where 7 G 



m \ 

1 + 2 a) 


and 0 decays like 


-( 1+2 a) 
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4.2 Consequences and remarks 

1. The speed of propagation for this two dimensional model cannot be purely exponen¬ 
tial. This also explains why trying to construct steady solutions to (4.3) does not lead 
anywhere. 

2. The influence of the road is felt on the expression of the solution to (5.1). Indeed, 
although there is no increase of matter on the road, the maximal growth rate - i.e. the 
same as in the ODE u = f(u ) - is chosen and the system behaves as if there were an 
effective growth term on the road. 

3. The effect of the mortality on the road is never felt in the growth exponent. Its 
only influence is that it divides the fundamental solution by a large factor. But, if one 
waits for a sufficiently long time, one will in the end observe propagation at exponential 
velocity. This suggests to study the asymptotics k —> +oo, and a transition of the type 
discovered in [26] for a —> 1. 

4. This raises the question of whether this is the correct asymptotics: sharp spreading 
rates are indeed not given, in general, by that of the linearized equation. The following 
simulation, taken from [ 5], investigates the issue. It solves the rescaled problem satisfied 
by v and u, defined on R x R + x R + , by 

v(x,y,t) — v(e lt t~ m x,y,t) and u(x,t) = u(e lt t~ m x,t). 

Here l = tww, and m ^ 0 the constant that we want to discover. 

l+2a 7 





Figure 4: Evolution of the density u with a = 0,5, for m — 0 (on the left), m = 2 {i+ 2 a) 
(in the center) and m = (on the right), at successive times t = 30,40, 50,..., 200 
with a color graduation from blue to red. 

The left side of Figure 4, that concerns m = 0, shows that the level sets move faster than 
ei+ 2a , whereas the right side, that concerns m = shows that the level sets move 

3 t _ Q 

slower than t x + 2a e 1 + 2a . The center of Figure 4 concerns the particular choice m = ^(\+ 2 a) > 
suggested by the upper bound of Theorem 5.1. On compact sets, the rescaled density u 
seems to converge to a function that does not move in time. 


17 



















5 The supersolution 

In this section, we are going to present the main lines of the proof of the 
Theorem 5.1 (Conlon | ]). Let (v(x,y,t),u(x,t)) solve 

d t v — Av = /'(0)h, i6l,j/>0,i>0 

d t u + (— d xx ) a u = —yu + v — ku, i6l,j/ = 0,i>0 

— d y v = yu — v, i6K,!/ = 0,t>0, 


(5.1) 


with (v(x, y, 0), u(x, 0)) = (0,Mo(ic)) and m 0 ^ 0 nonnegative and compactly supported. 
There exists a function R{t,x ) and constants 6 > 0, C > 0 such that 


1. we have, for large x: 


u(x,t) 


8a/isin(a7r)r(2a)r(3/2) e^A°)* 


n(k + /'(0)) 3 
2. and the function R(t,x) is estimated as 


t 3 / 2 |a;| 1+2 " 


^ R(t, x), 


0 ^ R(t,x) ^ C I e~ 5t + 


,/'(o )t 


of'{ o)t 


+ 


|/jr.|min(l+4a,3) | x |l+2a£ 


This is a computationally nontrivial result whose full proof will not be given here. 
However we will give the main steps. Let r 0 > 0 be the unique solution of 


rl = rl a + /'(0) + k. 


(5.2) 


will be denoted by r 0 . It is crucial to notice that r 0 > yjf'{ 0). 

We define (hi, Hi) = u ), the equation solved by (vi,ui) is 


a 


Vi 

Hi 


= -A 


Vl 

Hi 


—Avi 

(-a,)“Hi + yu x - Hi + ku\ + /'(0)Hi 


5.1 Fourier-Laplace transform 

The operator A is similar to the operator A defined in (A.l) with the constant k replaced 
by k + /'(0). Its domain is D(A) given by (A.2). It is a sectorial operator on X with 
angle /3j € (0, |). And so, we have 


Ui(x, t) 2i7T 1 y ’ 




0 

u 0 (x) 


e~ xt dX. 


(5.3) 


The computation of ( A — XI) 1 
gives 


u 0 


in the Fourier variables is an easy step that 


1 U ^ 


(A - XI) 


-i 


0 

u 0 


A 


p{\ lei) 


3 _ '\/~^+hF 




-a + KI 2 + i 
P( A,|ei) 


* Mq 


* Mo 




(5.4) 
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where P is given by 


F(A, |£|) ( — A + |£| “ + /i + k + f( 0 )) ^-A + |£| + 1^ — /r (5.5) 

The computation of this inverse Fourier transform requires the knowledge of the location 
of the zeroes of P ; one may prove that (see [21 ]): 

• if |£| < ro, for any A G C, P( A, |£|) does not vanish, 

• if |£| ^ ro, P(A, |£|) may vanish for some real values of A. 


5.2 First reduction 


The following preliminary lemma simplifies the expression of the inverse Fourier trans¬ 
form. Its poof is not given here, see [25]. 

Lemma 5.2. Let be defined in (5.2) and P be defined in (5.5). For r ^ 0, t > 1 and 
a constant (3 G (0, |) ; we set 


b (r,t) = — 


V-A + r 2 + l _ xt 


ITT 


- 0,0 


P(A,r) 


e dX, 


(5.6) 


where r 0j( g = © M+e 1/3 . Then, for all c G (0,1), we have 

1. For r G (0, cro) and t >1 : 

n —r 2 t r°° 


tpfat) = 


2 fie 


v 


e~ vt dv. 


7T J 0 |P(r 2 + u, r)\ 2 
2. There exists a universal constant C > 0 such that, for all r ^ cro a nd all t > 1 
Mr,t)\ < Ce^ r2a+k+f ' w - £o)t ^ \J |r 2 — (r 2 “ + k + f'(0) — e 0 )| + l) , 
where c 0 = ^“(l ~ c 2q ) > 0 . 


5.3 The main part {|£| < ro} 

Estimating the integral on {|£| < ro} is now a Polya type computation [ ], [ r ]. This is 

what we choose to develop here. 


Lemma 5.3. Define J(x,t ) as 

rcr 0 

J(x,t ) := / Ip A (r,t)e lxr dr. 

Jo 

where c e > 0 satisfy c 2 Vq cos(2c) > /'( 0). There is C > 0 universal such that, If |x| is 
large enough, we have 


J(x, t ) — 


4a;/xsin(a7r)r(2a)r(3/2) 

e f(o)t 

7T(k + /'( 0)) 3 

f 3 / 2 \x \ 1+2a 


^ Co ( g—c 2 rg cos(2e)t e -v^sin(e) _ 


X'- 
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Proof. From Lemma 5.2, we know that for r G (0, cro) and t > 1, 


tpA'pt) = 


2/ie 


—vt 


7T Jo | P(r 2 + v,r) 


du, 


where 

P(r 2 + zyr) = {-v -r 2 + r 2a + p + k + / / (0))(z v / ^ + 1) - p. 

We dehne, for (V, z) G M+ x C : 

Q(v, z ) = (- 1 / - z 2 + z 2a + k + /'(0)) 2 + z/(-z/ - z 2 + z 2a + /i + k + /'(0)) 2 , 
so that we have 

Q(z/, r) = |P(r 2 + z/, r)| 2 for (z/,r) G M + x [0, cr 0 ]. 

Tims, J becomes 


where 


J(x,t ) 


^ I'" 0 e- r2t e lxr j(r, t)dr, 
77 Jo 


j{r, t ) 



\fy 

Q(z/,r) 


e-^dz/. 


(5.7) 


To estimate J, there are three steps. The first one consists in rotating the integration 
line by a small angle e > 0. Then, we prove we can only keep values of r close to 0. 
Finally, we rotate the integration line up to |. 


Step 1 : From our knowledge of P, there exists a small angle e > 0 such that 

no) 

«*( 2<0 > ^ 

c r 0 

and 

for all v ^ 0 and z G {z G C | \z\ ^ cr 0 , arg(^) G [0,er]} : \Q{v,z)\ ^ cq. 

We want to rotate the integration line of £ in (5.7). For all t > 1, the function 

z ^ e~ zH e ixz j(z,t ) 


(5.8) 


(5.9) 


is holomorphic on the same set as Q, that is to say on {z G C | \z\ ^ cr 0 , arg(z) G [0, e]} 
if a G [|, 1) and on {z G C* | \z\ ^ cr 0 , arg(^) G [0, e]} if a G (0, |]. In this last case, 
we need to remove a neighborhood of zero when rotating the integration line. Choose 
5 G (0, cr 0 ). On the small arc 7 y e = {<5e* 0 , 0 G [0, e]}, we have for t > 1 



e ~ z tglXZ A ( 


j(z,t) dz^C e~ 5 ^ s mt e -xSsm(e) 6 d 0 : 


where C > 0 is a universal constant. The right hand side tends to 0 as 5 tends to 0. 
Thus, the Cauchy formula leads to 


J(x,t ) 


2/i 
7r 


(Ji( x ,t) - J 2 (x,t)), 


(5.10) 
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where 


rcr o 


J 1 (x 1 t)= I e- s2e2iSt+iE+ixseie j(se i£ ,t)ds, J 2 (x,t) = cr 0 i / e~*^ t+i0+ixcroeia jicroe*, t)d9. 
Jo Jo 


The term J 2 decays exponentially in time : 

\J 2 (x,t)\ < cr 0 j^ e -<?rlcos{26)t e -xcr 0 sin(0) | j ( C r 0 e iB , t) \ dO < £_ e -c 2 rlc OS {2e)t, ( 511 ) 


where C > 0 is a universal constant linked to Cq dehned in (5.9). 

Step 2 : We now treat J±. We cut it into two pieces in order to keep values of s close to 
0. Let us define, for x > (cro) -2 and t > 1 : 


r- 1 / 2 


J?(x,t) := I e~ s2e2iet e ixseis j{se i£ ,t)e i£ ds and 

Jo 


(5.12) 


r»cr 0 


J{(x,t):= I e~ s2e2iH e ixseis j{se ie ,t)e i£ ds, 
Jx - 1 / 2 


so that Ji(x,t) = J™(x,t) + J[{x,t). For x > (cr 0 ) 2 and t > 1, we have the estimate 


rcr 0 


\J[(x,t)\ ^ C / e -*ssin(e) e -s cos(2 

Jx - 1 / 2 


where C > 0 is once again linked to cq defined in (5.9). This implies that J{(x,t) decays 
exponentially in x and, taking C larger if necessary, for x > (cro) 2 and t > 1: 


\J[(x,t)\ ^ Ce~^ sia ( £ \ 


(5.13) 


Step 3 : We prove that J™{x,t) decays like x d+ 2Q ) f or large values of x. We turn 
the variable of integration into s = xs to get, for x > (cr 0 ) 2 and t > 1, 

rx 1/2 _ 2 r~ 

J™(x,t) = / e-& 2i ' t e aei 'j(sx- 1 e ie ,t)e ie —. 

Jo x 

Keeping in mind that we want an estimate for large values of x, we cut J] 71 as 


d / 2 




c 1 / 2 


J\ n (x,t) = / e lse j{sx~ 1 e l£ ,t)e l£ — + / (e 


-^r e 2ie t 


X 


1 )e isei£ j(sx~ 1 e i£ ,t)e i£ —. 

x 


The second term in the right hand side satisfies 


d / 2 


- 4 e 2 -t 




[e ~ - l)e“ e j(5x - V £ ,t)e 8£ — 

x 


<3 


^+00 


/y>3 I 

x J 0 


S 2 e - £sin(£) dx, 


(5.14) 
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where C > 0 is once again universal. We have to estimate 


r 1 / 2 


e iseis jisx-'e* ,t)e i£ — 


x 


for large values of x and t > 1. For all v G I 
Xo G (0,1) such that Q does not vanish in 
x~^ 2 < Xo, the function 


Q{v, 0 ) 7 ^ 0 . Consequenlty, there exists 
.+ x D Xq {{)). Thus, for all t > 1 and all 


z\-^e lz j(zx ,t) 


is holomorphic on {z G C | \z\ ^ x^ 2 } if a G [1/1) and on {z G C* | \z\ ^ a; 1 / 2 } if 
a G (0,1]. Let <5 G (0,1). On the small arc 7 s = {8e ld , 9 G [e, |]}, we have for t > 1 



e iz j(zx~\t)\dz ^ C 


e~ Ssin{e) 8d9. 


The right hand side tends to 0 as 8 tends to 0. For x > x 0 2 , we can rotate the integration 
line up to | and the Cauchy formula leads to 



e ise jCSx- 1 ^ 6 ,^— 
x 



d s 

e~ s j(isx L , t)i— 
x 


+ 



d/2 pie 


j(x 1/2 e ie ,t)ie 


iG 


(16 

x 1 / 2 ' 


The second term in the right hand side satisfies 


cV2ei9 j{x' 2 ^ ,t)ie id ^ 

0C ' 




(5.15) 


where C > 0 is a universal constant. It remains to estimate J / 1 dehned by 


d/2 


. ds 


J\ n (x,t) — / e S j(isx 1 ,t)i —, 

Jo % 


where 


j(isx ,t) = 


-e~ vt dv. 


J 0 Q{y,isx x ) 

Recall that we are interested in the real part of J{". A simple computation gives 


Q ( u, isx x ) 2 9fte 


Qiy, isx~ 


-)) = ^ isx ~' ) ^ m [o(^F0 


s 2 “ s 4q 

= 2 — sin(o; 7 r)((-z/ + s 2 x~ 2 + k + /'(0) + /i)(l + u) - /x) + — sin( 2 o; 7 r)(l + u). 


x 


x 


The integral under study is 


1/2 


Ke(JrOM)) = 


Am 


\ 


ds 


— y/Ue-«dv- 
Q(is, isx L ) J x 


With the dominated convergence theorem, we get 
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lim x 1+2Q 3f?e(J™(a:, t)) = 2 / e S s 2a sin(a;7r)h(z/, £)ds, 

ai-z+oo / n 


where h is defined by 

3 (—n + k + f(0) + /i)(l + z/) — /i 


h{y,t ) = 


This implies that 


iQM)r 


■\/ne ^dz/ ~ 


T( 3 / 2 ) 


z—h-oo (k + /'(0)) 3 f 3 / 2 


(5.i6) 

Finally with (5.13), (5.14), (5.15) and (5.16), we have the existence of a constant x\ > 
max(xo 2 , cr 0 , (cr 0 ) -2 ) such that, for all x > X\ and all t > 1 : 


jo,*) 


4asin(a7r)r(2Q;)r( 3 /2) 


(* + /'( o)) 3 


< C ^ e -^ sin(£) + ^ , 


with C > 0 universal. This estimate added to (5.10) and (5.11) leads to the existence of 
a constant C 2 > 0 such that, for x > x\ and t > 1, 

\J 2 (x,t)\ ^ C ^ e - c2r o cos(2e)i + e ^^sin( £ ) + _L^j , 


which proves Lemma 5.3. 


□ 


5.4 The remaining terms 

Since r 0 > \J /'(0), the integral on {|£| ^ r 0 } can be bounded from above by e _r ofi See 
[25] once again. 

6 An auxiliary subsolution for the transport equation 

Let us come back to our model (2.1). The following lemma provides a subsolution to 
a nonlinear transport equations with suitable exponential decay. In what follows, for 
A G M+, we set V\{x) := |x| -A . 

Lemma 6.1. Let g be a nonnegative function of class C°°(M), with g( 0) = 0, g'(0) > 0, 
and let a be a positive constant. Then, for all 0 < 7 ^ 7 := g'(fS)/a, the equation 

— r yxijj'(x) — g(if(x)), 16 R, ( 6 . 1 ) 

admits a subsolution 0^1 of class C 2 (M) and with the prescribed decay |(rp CT as |x| —y 00 . 
More precisely, there exist some positive constants {3, Ai, A 2 ,£ and D such that 

if |x| ^ A 2 , ~'yx(j)'-g((t>) < ~(3v a+£ , -0" A, Dv a+£ , (~d X3 ,)“0 ^ Dcf, (6.2) 

if |aj| G [A u A 2 ], - 7 X 0 ' - g((j>) < 0, (6.3) 

if |x| ^ A u 0 = 0(40. (6.4) 
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Proof. Let 5 G (0,1) be such that g is increasing on (0, <5) and, in addition, 


Vs G (0, S), g(s ) ^ </(0)s + (</'(0) - l)s 2 . (6.5) 


A first attempt to construct a subsolution satisfying the conditions stated in the lemma 
could be 


M x ) : = 


f V a (x) 

- Av a+e (x) 

if \x\ 

1 

A 

if x 

Ur 

A°+e 


with A, A\ > 0 to be chosen. Requiring that 0i G C 1 (M) yields 

A 1 = A l/e (l + £ -) l/S ■ 


With this choice, the function 0i is positive and nonincreasing on R + . But, if a ^ 1/2, it 
is not regular enough to yield an estimate of its fractional Laplacian. Consequently, we 
modify it to have a C 2 (M) function. This argument is, by the way, not so far from that 
of Silvestre in | | in the study of the regularity of solutions of integral equations. The 

function <fi is concave for Ai ^ |x| ^ A 3 , where 


A3 : — 


A l,£ 



£\ V e 

a J 



£ 

cr + 1 


l/e 

> Ai, 


otherwise it is convex. 



We fix a constant A 2 G {A\, A 3 ) and consider 

{ v a (x) - Av a+e (x) if \x\ ^ A 2 , 

x(M) if Ai < |x| < A 2 , 

Ja - = M A l) > o if |a?[ ^ A u 

where x G C 2 ([Ai, A 2 \) is nonnegative, nonincreasing, concave on (Ai,A 2 ) and chosen so 
that 0 G C 2 (M). First, we prove that 0 is a subsolution to (6.1) treating separately the 
different ranges of |x|. 
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A 


• |x| ^ A 2 ^ A 1/e . 

For A ^ 5~ £ / a , we have that 


0^0^ — 
Af 


A*+e 


< Ap < a _ ct/£ < 5, 


whence, by (6.5), g(0) ^ </(O)0 + (^(O) — 1)0 2 . This, together with 7 cr < </(0), yields 

-7 ~ g( 4 >) < 7CJ0-A7en ff+£ -5f'(O)0+(1-/(O))^ 2 

< (1-/(0))v 2 - A7£n CT+e . 


Since u 2 = n 2c7 , choosing £ < a, we have that the latter term is less than —f3v u+£ for 
given /3 > 0, provided A is large enough. This proves the first estimate in (6.2). 


• |x| G (A\, A 2 ). 

Using the concavity and monotonicity of y we see that the function r t —y 7 ry'(r) is 
nonincreasing in (Ai,A 2 ). Hence, by the monotonicity of g in (0,5), 

- 7 x 0 '(x) - g(<f>(x)) = - 7 xy'(|x|) - ^(y(|x|)) < - 7 A 2 x'{A 2 ) - s(y(A 2 )), 

which is less than or equal to — fdv a + £ {A 2 ) by the previous case. This proves (6.3). 

• [x| ^ A 1 . 

- 7 x 0 '(x) - g{<t>{x)) = -s(0(A)) ^ 0 . 

Then, we prove the remaining two estimates of ( 6 . 2 ) for |x| ^ A 2 . As for the first 
one, we have that 

-0"(x) ^ (a + e)(a + £ + 1 )Ax~ 2 v a+£ (x) ^ Dv a+e (x), 

for some D > 0. The last estimate, concerning (—d xx ) a (j), follows from the previous one, 
with a possibly smaller D. Indeed, as seen in [21], the function 0, being of class C 2 , 
radially symmetric and nonincreasing in |x|, satisfies (— d xx ) a fx)(p ^ ZT 0 (x) for some 
D' > 0. □ 


7 Lower bound close to the road 

Throughout this section, (v,u) is as in Theorem 2.2 and 777 * = /'(0)/(l + 2a). 
Namely, (v,u) is a solution to (2.1) starting from a nonnegative, compactly supported 
initial condition (u 0 ,Wo) ^ (0,0). The section is split into four subsections. In the 
first, short one, we estimate the solution of ( 2 . 1 ) at time to > 0 in strips of the form 
M. x [0, Y], In the second one, we study a reduced model - where we discard, from the 
full system ( 2 . 1 ), terms that seem to decay exponentially fast in time; we think that it 
will help in understanding the logic of the construction. The third subsection is devoted 
to the construction of a subsolution to the full model, which is then used in the following 
subsection to derive an auxiliary result to Theorems 2.2 and 2.3 part 2. 

Proposition 7.1. Under the assumptions of Theorem 2.2, for all 7 < 7 *, there holds: 

liminf inf u(x,t) > 0, (7.1) 

i^+oo Ix^e-V 

3A > 0, liminf inf v(x,y,t) > 0. (7.2) 

t —^-|-00 

O^y^A 
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7.1 Bounding from below the solution at positive time 

The next result provides a lower bound on the decay of solutions to the system with no 
reaction in the field (i.e., / = 0). 

Lemma 7.2. Let (p v ,p u ) be the solution to 


d t p v — A p v 

= 0, 

iGl,|/>0,l>0, 


d t p u + (-d xx ) a p u 

= -(/i + k)p u +p v , 

iGl,|/ = 0,l>0, 

(7.3) 

-d y p v 

= jip u — p v , 

iGl,|/ = 0,l>0, 



completed with the initial data p v (-, -,0) =0 and p u (-, 0) = uq. Then, for any constant 
Y > 0, there exists a > 0 such that 

V(x,y) Glx [0,y], p v {x,y, 2) ^ t^t + 2a > P u (x, 2) ^ “ ,i+ 2 a - ( 7 - 4 ) 

1 + \x\ 1 + \x\ 


Proof. Remark A.2 in the Appendix ensures that, for all time t > 0, the function p v is 
nonnegative on M. x M + . Thus the second equation of (7.3) gives 

Vi G 1, t > 0, d t p u (x , t) + (—d xx ) a p u (x, t) + (p, + k)p u (x, t ) ^ 0. 


Let us denote by p a the fundamental solution to the ID fractional Laplacian; it is well 
known that p Q decays like |x| _ ( 1+2 “) and a lower bound of p Q leads to the existence of a 
constant a\ > 0, depending on uq and a , such that for all and all t > 1: 


p u (x,t) ^ e *p a (x, t) ^ ai 


£ —i/2 ag — 


1 + lx 


1+2 a ‘ 


(7.5) 


Take xq large enough so that 


Wx G (x 0 - l,x 0 + 1), 


1+1 


l+2a 

1 + (1 + 

\x\) 1+2a 


1 

2 ' 


Set 


w(x , y, t ) := (1 + |z 0 | 1+2a )i/'(:co + x, y, t ). 


We have to prove that, for any Y > 0, there exists a > 0 independent of x 0 (large enough) 
such that w(0,y,2) ^ a for y G [0, Y]. The boundary condition satisfied by p v in (7.3), 
the estimate of p u in (7.5) and the choice of Xq give, for \x — Xq| < 1 and t > 1, 


—d y w(x,0,t)+w(x,0,t) = p(l + \x 0 \ 1+2a )p u (x,t) 

^ 2~ 1 a 1 t- 1 / 2a e~^ +k)t . 


If t G [1,3], this is larger than some constant 02 > 0. Thus, it easily follows from the 
uniform regularity of w with respect to y (which is the same as for p v ), that y t —y u>(0, y, 2) 
cannot be arbitrarily close to 0 at y = 0 without becoming negative at some y > 0, which 
is impossible. 

□ 
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7.2 A reduced model 

We are going to seek a stationary strict subsolution to the following system: 

- 7 ZdtV-dyyV = f(V), eeR,0 <y<L, 

-7 £d<:U = —yU + V — kU, £eR,y = 0, 

-d y V = yU - V, £ eR,y = 0. 


(7.6) 


for a large constant L. We want the subsolution to have the algebraic decay |£| ^ 1+2a ^ 
for large values of |£|- Taking 

L> tt (/'(0) - (1 + 2 a) 7)“ 1/2 = 7 T [(1 + 2a)( 7 * - 7)]~ 1/2 , (7.7) 

we see that 7 < (/'(0) — n 2 /L 2 )/(l + 2a) and hence Lemma 6.1 applies with 

g(s) := f(s) - (^) s, a 1 + 2a, (7.8) 

providing us with a function 0 satisfying (6.2)-(6.4) and decaying as \x\~^ 1+2a \ Define 

0(0 sin (J^y + /i) if 0 ^ y < L (l - J) 

0 ify^L(l-f) 


H00 = 


mo = cm , 


where 0 < h < n and C > 0 will be suitably chosen in such a way that (V_,U) is a 
subsolution to (7.6). We treat separately the three equations of the system. 

• The first equation. 

The nontrivial case isO<?/<L(l — -), and there holds 


- d„V - f(V) = f-7 w + (£) 2 A sin Qs + ft) - f(V) 

= - 7 ^ 0 'sin (Oy + fi) - g [(p sin (Oy + 0 )) 

< (—7^0' - 1/(0)) sin (Oy + h) , 

where the last inequality holds because s 1 —> g(s)/s is decreasing. Then, using the 
properties (6.2)-(6.4), we derive the existence of /3,e > 0 such that 

sin (J^y + h^j . (7.9) 


-7 mm-dyyV-nv)^- 

The second equation. 

Since —7£0' 0 g(0) =7 0(0)0, we have that 


0 


1 + 10 


1 + 2 a+e 


- 7 iU! + + k)U - HO 0) = C[- 7 ^ 0 ' + (y + fc)0] - 0sin(/i) 

< [(7(0(0) + y + A:) — sin(/i)]0. 

Consequently, for C = C(/i) small enough, 

- 7 OW + H + fc)C7 - HO °) < -00, 


(7.10) 


for a possibly smaller (3. 
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The third equation. 
There holds 


- d y V(£, 0) - nU{£) + He, 0) = (-^ cos (h) - Cp + sin {hj) 0(0 < 0, (7.11) 


up to choosing h < arctan(7r/L) and then C = C(h ) small enough. 

7.3 Subsolution to the full model (2.1) 

The subsolution (v, u) that we are going to construct is a modification of the pair (H U_) 
defined in the previous subsection. We know from Lemma 6.1 that (V_,U_) decays as 
|£| -1-2 ", an d H (6-2), that there exists D > 0 such that 




Lemma 7.3. For B > 0 small enough, the couple (v, u) defined by 
v(x,y,t) — V(xb(t),y) and u(x,t) — U_(xb(t)) : 


(7.13) 


with b(t) = Be ^, is a subsolution to (2.1). 

Proof. Let us call 

£i{v) = d t v - An - f(v), £ 2 (v,u) = d t u + ( -d xx ) a u + (n + k)u - 7 0 v. 

• In the field. 


The nontrivial case is 0 < y < L( 1 — -), where, by properties (7.9) and (7.12), we get 
A(l) = - 7 xb(t)d^V(xb(t),y) - b 2 {t)d^V_(xb(t),y) - d yy V_(xb(t),y) - f (V(xb(t ), y)) 



Then, since \b(t)\ ^ B, this term is negative for B small enough. 


• On the road. 

Using (7.10) and the second estimate in (7.12) we derive 


C 2 (y, u) = — 7 xb{t)Uf{xb{t)) + b 2a (t)(—d^) a U_(xb(t)) + (y + k)U_{xb{t )) 
-V(xb(t), 0 ) 

< (~j3 + B 2a D)U(xb(t )) 


which is negative for B small enough once again. 

• The boundary condition is an immediate consequence of (7.11). 
This shows that (y,u) is a subsolution to (2.1). 


□ 
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7.4 Conclusion by comparison with the subsolution 

Proof of Proposition 7.1. We derive the lower bound close to the road by fitting the 
subsolution (v, u) provided by Lemma 7.3 below (v, u ) at time t — 2. To do this, we 
make use of the pair (p v ,p u ) from Lemma 7.2. Recall that, by construction, 

C C 

v(x,y,2) C i+ 2 a l[o,£](y), u(x,2) < , 1+2a - 

1 + \x\ 1 + |x| 

for some L, C > 0. Since the nonlinearity / is nonnegative, the comparison principle 
of Theorem A.2 entails that (v,u) is greater than (p v ,p u ). Applying Lemma 7.2 with 
Y = L, we infer the existence of a constant a > 0 such that for all i 6 l and y E [0, L ], 

v(x,y,2) ^p v (x,y,2) > -- ° 1+2a , and u(x, 2) ^ p u (x, 2) ^ Q 1+ , a . 

1 + |x| 1 + \x\ 

We eventually infer that, at time t — 2, (v,u) is greater than £o(h, Ah provided £ 0 > 0 is 
small enough. Notice that since s (->• is decreasing, for e 0 E (0,1) the couple s 0 (v,u) 
is still a subsolution to the problem (2.1). Therefore, choosing £o £ (0,1) sufficiently 
small, we can apply the comparison principle and obtain 

W(x,y) 6 ix M + , t ^ 2 , v(x,y,t ) ^ £ov(x,y,t), u(x,t ) ^ £q u(x,t). 

Finally, we know from Lemma 7.3 that 

v(x, y, t) = V^fxBe^, y ), u(x, t ) = UfxBe ' 7 *), 

with B > 0 , V positive in some strip M x [0, A] and U_ positive. The proof of Proposi¬ 
tion 7.1 is thereby achieved. □ 

8 Lower bound in the field 

This section is dedicated to the proof of following weaker version of Theorem 2.3 part 2: 
Proposition 8.1. Under the assumptions of Theorem 2.2, there holds 

\/9 E (0, 7 r), 0 < c < ck/suYO ), liminf inf v(r cos(60, rsin( 6 ) ),t) > 0. 

t —>+oo O^r^ct 

Proof. As said in Section 2 , the invasion on the road is exponential in time, whereas it 
cannot be more than linear in the held. Therefore a good model for it is the Dirichlet 
problem 

d t v - dyyV = f(v) for y > 0, v( 0 , t) = 1 . 

Extend v as a function of two spatial variables by v(x,y,t) := v(y,t ): this gives the 
desired propagation. In order to make this consideration rigorous, we first exploit the 
lower bound close to the road given by Proposition 7.1, next we use standard arguments 
for the spreading in the held orthogonally to the road. 

Let (v,u) be as in Theorem 2.2. By (7.2), for any 7 E ( 0 , 7 *), there exists 6, to, A > 0 
such that 

Vi ^ to, |x| ^ e 7t , 0 ^ y ^ A, v(x,y,t ) ^ 8. (8.1) 

We now forget the road and replace it with sudden death of the population, namely, with 
the Dirichlet boundary condition. We derive the following 
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Lemma 8.2. Let w be the solution to the problem 


J d s w — Aw = f(w), iGl, y > 0, s > 0, 

0, s) — 0, x G I, s > 0, 

starting from a nonnegative, compactly supported initial datum Wo ^ 0. Then 
VO < c < ck, lim w(x, y + cs, s ) = 1, Vc > ck, lim w(x, y + cs, s ) = 0, 

S—>-+00 5—>-+00 


( 8 . 2 ) 


(8.3) 


locally uniformly in (x, y) el 2 . Moreover, for any K d M, p > 0 and 0 < c < ck, 


inf v(x, y, s) > 0. 

s^p, x£K 
p^y^.cs 


(8.4) 


It is also possible to show that w converges locally uniformly to the unique positive 
bounded solution W to the ODE —W"(y) = f(W(y)) for y > 0, such that W(0) = 0. 
Let us postpone the proof of this lemma until the end of this section and continue with 
the proof of the proposition. We actually derive a stronger result: the uniform lower 
bound in rectangles expanding with any speed less than Ck in the vertical direction, and 
exponentially fast in the horizontal one (cf. Figure 4). 


y 


csin(0)£ 



-e 7 £t 


e 7 £t 


x 


Figure 4: Expanding rectangle 


Fix c G (0, Ck ), 0 G (0,7r) and £ G (0,1). Let t\ ^ t 0 /e be such that 


Vt > t u e 1£t ^ c ' C ° S (h!' t + 1. 

sm(0) 

It follows from (8.1) that, for all t^t\, there holds 

V|xq| ^ e 1£t — 1, |x| <1, j/G [0, A], v(xq + x,y,et) ^ 6. 


(8.5) 


Consider w the solution to (8.3) with initial datum w(x,y, 0) = 51(-i,i)x(o,a)(25 y)- Since 
v is a supersolution to (8.3), the comparison principle yields 

Vt ^ t\, |xo| ^ e 1£t — 1, igK, y ^ 0, v(xq + x , y , t ) ^ w(x, y, (1 — e)t). 

Applying the estimate (8.4) from Lemma 8.2 with x = 0 and s = (1 — e)t, we derive 


inf v(x o ,y,t)>0, 

p^.y^(l-e)ct 

for any p > 0. This estimate actually holds true up to p = 0 thanks to (8.5). Notice 
that, for 0 ^ r ^ sin c ^ (1 — e)t, we have |rcos(0)| ^ e 1£t — 1 by the choice of t\, and 
r sin(0) ^ c(l — e)t. The proof of Proposition 8.1 is thereby complete owing to the 
arbitrariness of c G (0, Ck) and £ G (0,1). □ 
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Proof of Lemma 8.2. The second limit in (8.3) follows immediately from the spreading 
result of [3, 40], because solutions to (8.2) are subsolutions of the problem in the whole 
plane. For the same reason we know that lim sup s _^ +0O w(x, y, s) ^ 1 uniformly in (x, y) G 
M. x M + . Let us deal with the fist limit in (8.3). Fix c G (0, Ck) and cast the problem in 
the frame moving vertically with speed c. That is, consider the problem for w(x,y, s) := 
w(x, y + cs, s): 

| d s w - Aw - cdyiv = f(w), iGM, y > -cs, s > 0, , g ^ 

1 w(x, —cs, s) = 0, x G !, s > 0. 

We need to prove that lim inf s _>. +00 wfx, y, s) ^ 1 locally uniformly in (x,y) G M x M + . 
Let A c (r) be the principal eigenvalue of the operator —A — cd y in the two-dimensional ball 
B r , under Dirichlet boundary condition, and ip c be the associated positive eigenfunction. 
This operator can be reduced to a self-adjoint one by multiplying the functions on which 
it acts by exp(— (c/2)y). This reveals that A c (r) — c 2 /4 = Ao(c), the principal eigenvalue 
of —A in B r . Hence 

2 2 

lim A(r) = lim A 0 (r) + ~r = ~r < f'( 0). 

There is then r > 0 such that /(s) ^ A(r)s for s > 0 small enough, and we can therefore 
normalize the principal eigenfunction ip c in such a way that, for all n G [0,1], Kip c is a 
stationary subsolution to (8.6) for (x, y) G B r and s > r/c. Moreover, for given So > r/c, 
up to reducing k if need be, K,tp c lies below the function if at a time s = s 0 , the latter 
being positive by the parabolic strong maximum principle. Let w be the solution to 
(8.6) emerging at time s = Sq from the datum mp c extended by 0 outside B r , which is a 
generalized subsolution. It follows that w is increasing in s and converges, as s —» +oo, 
locally uniformly to a positive bounded solution W of 


—AW - cd y W = f(W) in M 2 . 

Then, from one hand, W = 1 by the Liouville-type result of [ 1], Proposition 1.14. From 
the other, the comparison principle yields w ^ w for s ^ so- AVe eventually derive that 
w —> 1 as s —» +oo. 

We finally turn to (8.4). From the above arguments we know that, for R > 0 large 
enough, the principal eigenfunction <p 0 of — A in Br satisfies —A<p 0 ^ f(<p o)- By B ie fist 
limit in (8.3) we have that, up to renormalizing <p 0 if need be, for given c G (0, ck), there 
exists so > R/c such that 

Vs ^ s 0 , {x, y) G B r , w(x, y + cs, s) ^ <p 0 (x, y). 

Hence, by the comparison principle, 

Vs ^ s 0 , t ^ 0, (x, y ) G Br, w(x, y + cs,s + t ) ^ ip 0 (x, y ). 

Consider r ^ So and cso ^ r) ^ cr; applying the previous inequality with s = rj/c (^ so), 
t = t — rj/c (^ 0), y — 0, we infer that 

inf w(x, ii, t ) ^ inf <p 0 (^, 0) > 0. 

T^ 0 , \x\^R/2 \x\^R/2 
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By the parabolic strong maximum principle, in the above infimum, r A s () can be replaced 
by r ^ p, for any given p > 0. We can further get a positive lower bound for w on the 
set r ^ p, |x| ^ Rj 2, p ^ r) ^ csq by comparison with the function ip 0 (x,y — R — p/2), 
suitably normalized. Namely, it holds true that 

inf w(x,ij,r)> 0, 

from which (8.4) eventually follows by a covering argument, owing to the invariance of 
the problem by x-translations. □ 

9 Convergence to the steady state in the invasion set 

In this section, we put together the previous results and derive Theorems 2.2 and 2.3. 

Proof of Theorems 2.2 and 2.3 part 1. We start with Theorem 2.2. Let (v, u) be as there. 
The pair (v,u) given by Theorem 5.1 is nonnegative by the comparison principle given 
in the Appendix, and thus it is a supersolution to (2.1) because / is concave. The same 
holds true for K(y,u), for any K > 0. We choose K large enough so that, at time t — 1, 
K(v, u) is above the compactly supported initial datum of (v, u ). This is possible because 
v, u are strictly positive at any time t > 0, as is readily seen by applying the parabolic 
strong maximum principle and Hopf’s lemma to derive the positivity of v, and next using 
the equation for u. Therefore, by comparison, (v,u) lies below K(v(-, ■ + 1 ),«(-, • + 1)) 
for all times and thus Theorem 2.2 part 1 follows from Theorem 5.1. 

We now turn to Theorem 2.3. Let / be a concave function vanishing at 0 and at 
IMloo + 1 and such that 

7 >f in (0, IML + !)» 7(0) = /'(0). 

Consider the solution to 

d f v - d yy v = J(v), y e M, t > 0, 

starting from any positive bounded initial datum. We know from Aronson and Weinberger 
[ ] (or even [ ]) that v spreads at speed ck = 2a J /'(0), that is, 

Vc < ck, lim inf v(y,t ) = ||u|| +1, Vc > ck, lim sup v(y,t ) = 0. (9.1) 

t ~>+°° t^ + OO 

From the first property (recalling that no is compactly supported) we deduce in particular 
the existence of to > 0 such that 

VyeR, v(y,t 0 ) ^ v 0 (y), Vt ^ t 0 , v(0, t) ^ IMU ■ 

It follows that v(y,to + t) > v(x,y,t ) at t — 0 and also for t ^ 0, x G M, y = 0. 
Namely, v and v(y, t 0 + 1) are respectively a solution and a supersolution to d t — A = / 
in ( x,y,t ) € M x M + x M + , which are ordered at t — 0 and at y — 0. The comparison 
principle eventually yields 

V(x, y) G M x M + , t ^ 0, v(x, y, t ) ^ v(y, t 0 + t). 

Statement 1 of Theorem 2.3 then follows from the second condition in (9.1). □ 


32 



Proof of Theorems 2.2 and 2.3 part 2. Fix 7 G (0,7*) and let (x T ) T>0 in M and (y T ) T >0 
in M + be such that 

\x r | ^ e 7T , (y T )r> 0 is bounded. 

It follows from the estimates in the Appendix that, as r —» +00, the family 

(v(- + x T , •, • + r), u(- + x T , • + r)) T> 0 

converges (up to subsequences) locally uniformly to a bounded solution (fu, u) of (2.1) 
for all t 6 1. We claim that (v, u) coincides with ( V S ,U S ), the unique positive bounded 
stationary solution. This would yield Theorem 2.2 part 2. From (7.1) in Proposition 7.1 
(applied with a slightly larger 7) we deduce that 

m := inf u(x, t ) > 0. 

(ai,t)GlR 2 

Then, on the one hand, by Theorem 2.1, the solutions (Vi, Uf) and (V 2 , U 2 ) to (2.1), with 
initial datum (0,m) and (HnHoo, ||u||oo) respectively, tend locally uniformly to ( V S ,U S ) as 
t — > +00. On the other, by comparison with (n(-, •, •— n), u(-, •—n)), for all ( x , y) G lxR + 
and n G N, we have: 

Vi ^ 0, (W(^,y,i),^i(a;,i)) ^ (v(x, y, t - n), u(x, t - n)) ^ (V 2 (x,y,t),U 2 (x,t)), 

whence, calling s = t — n and letting n —> +00, we derive 

Vs G M, (v(x, y, s),u(x, s)) = (V s (y), U s ), 

that was our claim. 

It remains to prove Theorem 2.3 part 1. Fix 6 G (0,7r), c G (0, Ck) and consider a 
family (r T ) T> o such that 0 ^ r T G^ cr. We need to show that 

v(r T cos(9),r T sin(6),T) — V s {r T sin(6 1 )) —>• 0 as rG +00. 

This has been done above for sequences of r for which r T is bounded, because x T := 
r r cos(0) ^ e 7T , for any 7 > 0 and r large enough, and y T := r T sin(6 ) ) is bounded. 
Consider the case r T —> +cxd as r —» +cxd. The sequence of translations (v(- + r T cos (9), • + 
r T sin(0), - + r) converges (up to subsequences) locally uniformly as r —» +00 to a bounded 
function v satisfying 

d t v — Av = f(v ), (x, y) G M 2 , t G M. 

Furthermore, applying Proposition 8.1 with a slightly larger c and values of 9 in a neigh¬ 
bourhood of the fixed one, we infer that inf v > 0. Thus, by comparison with solutions 
of the ODE V' = f(V), one readily gets v = 1, which concludes the proof of Theorem 
2.3 part 2 because IA(+oo) = 1. □ 
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Appendix: Cauchy Problem and comparison principle 

We choose to prove existence, uniqueness and regularity to the Cauchy Problem for 
(2.1) by the theory of sectorial operators and abstract theory of semilinear equations, as 
exposed in Henry [37]. Comparison is then proved by standard integration by parts. We 
point out that this is not the only way, we could also use viscosity solutions theory. This 
would not, however, yield a significantly shorter study. 


A.l Existence, uniqueness, regularity 

We work in the Hilbert space X = {(v,u) G L 2 (R x M + ) x L 2 (M)}. This framework will 
be sufficient for what we wish to do. In what follows, 70 and 71 denote the usual trace 
and exterior normal trace operators. The operator A is defined by 



—Av 

(-d xx ) a u + nu - 70U + ku 


(A.l) 


The domain of A is 


D(A) = {(v,u) G H 2 (R x M+) x H 2a (R ) | 7m = pm- j 0 v} C X. (A.2) 

Notice that the first component v of an element of A is a continuous function by the 
embedding result, thus the trace 70 u is simply the value of v at y = 0. From [ 5], the 
operator A is sectorial in X. Cast the problem (2.1) in the form 


W t + AW = F(W), where W 



f(v) 

0 


J 


where / is C°°(M) and extended in a smooth fashion so that /(0) = /(1) = 0, and / = 0 
outside (—1, 2). From Theorem 3.3.3 and Corollary 3.3.5 of [ 7], there is a unique global 
solution (v,u) to (2.1) starting from (w 0 ,Wo) £ X such that 1 1—> (v,u)(-,t) belongs to 

C((0, +00), H 2 (R x M + ) x H 2 a (R)) n C^+oo), L 2 (R x R+) x L 2 (M)). (A.3) 

We are going to prove the 

Theorem A.l. The solution (v,u) of (2.1) starting from (0,Wo), f or a continuous, non 
negative and compactly supported function uq ^ 0 , satisfies (■ v,u ) G C°°(M x M + x M^_) x 
C°°(M x R* + ). 


Proof. This regularity result will be obtained by induction, proving the existence of a 
constant 6 G (0,1) such that, for all T > 0, e G (0, T) and n G N : 

v e C n+5 ’^(R x R + x (£ 2n ,T]) and u G C n+5 ’^(M x (e 2n , T}), (A.4) 


2 n 

where £ 2 n = / 2 _J £ —> e. 

‘ J n^+oo 

3 =1 

We will use the following results : 

Result 1 : Elementary Sobolev embeddings. We have : 
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• H 2 (R x K+) C C A (M x M + ), for all A G (0,1), 

• for a G Q, 1) : H 2a (R) C C 2 “-5(M). 

Result 2 : Let 0 < to < T, l > 0, and consider two functions g G C 1 ’* (M. x M + x [to, T ]) 
and U\ G (R x [t 0 ,T]). Then, from Theorem 4.5.3 in [ 1], the solution V\ of 

d t v x - Avi = g, x G R, y > 0, t > t 0 , 

—d y v i + i?i = gu i , iGl,|/ = 0,i>to, 

starting from Vi(-, t 0 ) G C i+2 (M x M+), for which the following compatibility conditions 
of order rri\ := [^J hold : 


d[ m \-d y v i + ui)(-, -,t 0 ) = /n9 t (m) ni(-,t 0 ), for all m = 0,... ,m 1 , 
satisfies G C i+2, 5 +1 (M x M + x [t 0 ,T]). 

i 

With these two results, we can prove (A.4). Let T > 0, e G (0, T) and £ t = 2~ J e. 

3 = 1 

In view of Result l,a = l/4isa special value. Thus we break the study into two cases: 
a ^ 1/4 and a > 1/4. We detail the latter case, and explain the needed modifications 
for the former one. 

Case 1. a > 1/4. 

• Case n — 0 : Since (v,u) belongs to the space in (A.3), Result 1 and Lemma 3.3.2 of 

[ | yield a constant <5 > 0 such that (v,u) G x M + x (0, T]) x x (0,T]). 

• Case n = 1 : We first prove that (d t v,d t u) G C((Ai, T]), H 2 (R x M + ) x H 2a (R )) D 
C 1 ((ei, T]), L 2 (M x R + ) x L 2 (M)). It is sufficient to prove that (d t v,d t u) is solution to 

d t w + Aw = F\ (w, t ), t > e i, (A.5) 

starting from (dtv(-,-,£i),d t u(-,£i)) G L 2 (M x M + ) x L 2 (M), where F\ is defined on 
A" x M + by 


F\ 



w l 

w 2 



0 I 


(A.6) 


As is usual, we can not directly differentiate equation (2.1) with respect to time, that 
is why we consider, for h > 0, the functions Vh and Uh defined on M x M + x M + by 


+ h) u(-,t + h)-u(-,t) 

v h (-,-,t>-= - t - and ■■= - - -. 


For any h > 0, ( Vh,Uh ) is in D(A) and satisfies 

V h \ + A ( V h \ = ( «v.at))-/Mv,<)) 

u h 1 V Uh ) \ o 


, t > 0. 


(A.7) 
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Once again from Lemma 3.3.2. of [37], v(x, y, •) is Holder continuous in time, uniformly 
in (x, y) E Rx K + , which implies that F\ satishes the assumptions of Theorem 3.3.3 and 
Corollary 3.3.5 of [ |. Thus, from [52], we can pass to the limit as h tends to 0 in (A.7) 

to get that (d t v,d t u) is the solution to (A.5), starting from (d t v(-,-,£i),d t u(-,£i)) G 
L 2 (R x M + ) x L 2 (R). Thus, we conclude 

(d t v, d t u) G C((e i, T],H 2 (R x R + ) x H 2a (R)) nC 1 ]^, T], L 2 (R x R+) x L 2 (R)). (A.8) 

We now study, for all t G [£2,7"], the couple (d x v(-, -, t), d x u(-, t)). We first prove that, 
for all t G (0, T\, d x u(-,t ) exists and 

d x u(-,t) G L 2 (R). 

From (A.3) and (A.8), we know that, for all t G [£ 2 ,T] : 

u(-, t) G H 2a (R ), d t u(-, t) G H 2a (R) and v(-, -, f) G ih 2 (M x M+). 

Applying the operator (— d xx ) a to the equation 

d t u(x,t ) + (-d xx ) Q u(x,t) = -{y + k)u(x,t) +v(x,0,t), 
we have for all f G [£2,7"] : 


(- d xx ) 2a u(-,t ) = -(/x + fc)(-9 !E!E ) a M(-,t) + (-9 ;K;r ) O! n(-,0,f) - (-d xx ) a d t u(-,t). 


This proves that for all t G [e2,T], u(-,t) G 7/ 4 "(M) C As done in the case 

n — 1, it is sufficient to prove that (d x v,d x u) is the solution to (A.5), starting from 
( d x v( -, •, £ 2 )j d x u(-, £ 2 )) £ L 2 (M x M + ) x L 2 (M), where is defoied in (A.6). 

Once again, we can not directly differentiate equation (2.1) with respect to x, that is 
why we consider, for h > 0, the functions Vh and Uh defined on M x M + x M + by : 


Vh{x,-r) 


v(x + h,-,-) -v(x,-,-) 
h 


and Uh(x, •) 


u(x + h, •) — u(x, •) 
h 


Passing to the limit as h tends to 0 in the problem solved by ( Vh,Uh ), [52] gives that 
(d x v,d x u) is solution to (A.5) with ( d x v (■, •, £2), d x u(-, £ 2 )) as initial datum, and conse¬ 
quently 

(d x v, d x u) G C((e 2 , T],H 2 (R x M+) x H 2q (R)) n ^((£2, T], L 2 (R x M+) x L 2 (R)). 


• Case n = 2 : To get the regularity of n, we apply Result 2 with 

wr =« G C 1+ ^(M x [£ 3 ,T]), g = f{v) G C 1+<5 ’^(M xl + x [£ 3 ,T]), 

and initial condition n(-,-,£ 3 ), to get that v G C 2+5,1+ ^(M x M+ x [£ 4 ,T]). It remains 
to prove the regularity on u, more precisely that d t u G C 5, 2(M x (£ 5 ,T]) and d xx u G 
C <5, 5(M x (£5, T}). This is done as in the case n = 1, applying several times the operator 
{—dxx) 1 to the equation for u, 7 being any positive number < a. 
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Iterating, we get (A.4). 

Case 2. a ^ 1/4. It is enough to show that there is $ > 0 such that (v,u) G C 5, 2(M x 
M + x (0, T]) x C s ’% (M. x (0, T]). From that, the proof of Case 1 applies. To prove that, one 
applies alternatively Theorem 3.3.3 and Corollary 3.3.5 of [ ] to get, inductively, that 

• u G C((0, T), H ka (M.)) and v G C((0, T), H 3 / 2+ka (R 2 + )) 

• d t u G C s ((0,T),H^ Q (R)) and d t v G tf3/2+(fc-i)a( M 2j) 

• and then, by Lemma 3.3.2 of [37], that u G C 5 ((0, T), iC fc_1 ) a (M)), 

• and, finally, that v G C 5 ((0, T), if3/2+(k-i)«( R 2 

We are back to Case 1 as soon as ka > 1/4. □ 

Once we know that the solution (v,u) to (2.1) is regular in space and time, we can 
remove the trace operators and the Cauchy Problem is thought of in the classical sense. 

Remark A.2. In the particular case of ^-independent solutions, a similar proof as the one 
done implies the unique solvability in C((0, +oo) , H 2 (R + ) x M) fl C((0, +oo) , L 2 (R + ) x 
M) for (2.1) with an ^-independent initial condition in L 2 (M + ) x M. 


A.2 Comparison principle 

It will follow from standard arguments; however we give some details: it is a crucial tool 
in the whole. Recall first that, from a straightforward computation, we have 

/ (—A ) a h(x)h + (x)dx ^ 0, 

Jr 

for all h in H 2a (R), 


Theorem A.3. Let(v\,ui), (i’2,1*2) be two couples in C((0, +00 ), (R 1 (M x M + ) flLip(M x 
M + )) x H 2a (M.)) fl C 1 ((0, +cx)), L 2 (M x M + ) x L 2 (M)) that satisfy 


and 



- 

Aui 



+ ( _ <9x3 

:)"«! 




7Ri 



dtV2 ~ 

Au 2 



+ (~d xx 

:) a «2 




7l' y 2 



f(v 1), 

-/aui + 7 0 fi - kui, 

hW - 7o<h, 

f(v 2 ), 

—p,U2 + 7cd ; 2 — ku2, 
IIU 2 - 7o'«2, 


iGK,!/>0,1>0, 
iGK,!/ = 0,1>0, 
iGR,j/ = 0,t>0, 

x G M, y > 0, t > 0, 
x G M, y = 0, t > 0, 
x G M, y = 0, t > 0. 


If for almost all ( x , y) G M x M + v±(x, y, 0) ^ i’ 2 (x, y , 0) and ui(x, 0) ^ U 2 {x, 0), then for 
all ( x , i/,1)gKx M+ x M +; we have vi(x, y , t ) ^ ^(x, y , t ) and ui(x, t ) ^ U 2 (x, t). 


Proof. Let / > 0 be a constant greater than the Lipschitz constant of /. Set 


(v 3 (x,y,t),u 3 (x,t)) := (ui(x, j/,t), wi(x,t))e 14 - (n 2 (x, y, t), u 2 (x, t))e lt , (A.9) 
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we have 


d t v 3 - Av 3 
dtu 3 + {—d xx ) a u 3 
l\v 3 




e~ lt f(v i) - e~ lt f[y 2 ) - Iv 3 , 
-/i« 3 + 7 0 f 3 - ku 3 - lu 3 , 
iiu 3 - 7 0 u 3 , 


iei,po,t>o, 

i 6 l,|/ = 0 ,t> 0 , 
i6t,|/ = 0,i>0. 


(A.10) 


Almost everywhere in R x R + and R we have u 3 (-, •, 0) ^ 0 and w 3 (-, 0) ^ 0. Multiply the 
first equation of (A. 10) by and integrate over R x R + to get 


d t v 3 vjdxd/y — // An 3 n 3 dxdy ^ // (Z |u 3 | «- lv 3 )v% dxdy — 0. 

f x£M.,y >0 «/ J 0 J J x€M.,y >0 

(A.H) 

Since rq and n 2 belong to C((0, +oo), //’ 1 (R x R + )) D C 1 ((0, +oo), L 2 (R x R + )), we have 
classically 

d t v 3 v$dxdy = ( j J \v£\ 2 dxdy) . (A.12) 

lxeM,y>0 ^ \JJxeR,y >0 J 


Using the third equation of (A. 10) and the fact that 70U3" ^ 0, we have 


A v 3 v^dxdy^— // |Vu^| dxdy + y / u^ov^dx— / |7ou^| ofe. 

a?GM,2/>0 */ «/ x£M.,y>0 JM. Jm. 

Inserting this last inequality and (A. 12) in (A. 11), we get 


1 d 

2 dt 


\v^\ 2 dxdy) / / I Vuj"| 2 dxdy + [i I u^ov^dx. 


' xE.^,y >0 / J J xGM,y>0 «/R 

Working similarly with we get 

L uhavtdx - 

The continuity of the trace operator gives a constant C tr > 0 such that 

Y2 


11 70^3” 


1 2 < C 2 llbi+ll 2 ||Vu + l | 2 

Il 2 (R) ^ tr U| l 3 |Il2(RxR + ) ^ 11 V 6 3 |Il 2( RxR+ ) 


Add (A. 13) and (A. 14), and use the continuity of the trace operator: 


|Il 2 (RxR + ) 


(A.13) 


(A.14) 


1 d /11 _j_ I j 2 11 112 \ /11 _|_ I j 2 .. _ 

2dt V^ 3 ll- L2 ( K ) l^ 3 lli 2 (®xR + ) J ^ ^ \W U 3 IIl 2( R ) + \\ V 3 II 7-,2n»vR ) I ’ 


^3 IIl 2 ( R ) T ||^ 3 || L 2 ( RxR+ y 

where C > 0 is a universal constant. Since w 3 (-,0) = 0 and v^(-,-,0) = 0 almost 
everywhere, we have for all t ^ 0, u^(-,t) = 0 and = 0 almost everywhere, 

which concludes the proof. □ 


The above comparison principle is stated for classical solutions whose initial condition 
belongs to X = L 2 (R x R + ) x L 2 (R). However, it is necessary for later purposes to have 
a similar result for solutions starting from ^-independent initial data, which therefore do 
not belong to X. 
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Theorem A.4. Let (vi,u\) be as in Theorem A.2, and ( 112 , 112 ) € C((0, + 00 ), H 2 (R + ) x 
M) fl C 1 ((0,+ 00 ), L 2 (M + ) x M) be such that 


f d t V 2 ~ dyyV 2 


f(v 2 ), 

i6l,i/>0,l>0, 

l u 2 


—/i «2 + 7 o '«2 — ku 2 , 

i6l,t/ = 0,l>0, 

{ llV 2 


tm 2 - 7 0 u 2 , 

i6l,t/ = 0,l>0. 

Ifvi(-, -,0) ^ V 2 (y,0) and u\( 

:-,o) 

^ ^ 2 ( 0 ) , then ni(-, 

•, t ) ^ v 2 (-, t ) and ui(-, t ) ^ u 2 (t) 


The proof is a straightforward adaptation of that of Theorem A.2. indeed, the points 
where the integrability is required involve the positive parts of the couple (vs,Us) defined 
by (A.9), which do belong to the desired spaces. For the same reason, one can also 
handle the case where both the sub and the supersolution are x-independent and the 
subsolution belongs to C((0, + 00 ), H 2 (R + ) x M) flC 1 ((0, + 00 ), L 2 (M + ) x M), which is even 
simpler because one is reduced to a problem in one less spatial dimension. 
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